Search for: All records

Abstract GW231123, the most massive binary black hole (BBH) merger detected by LIGO–Virgo–KAGRA, highlights the need to understand the origins of massive, high-spin stellar black holes (BHs). Dense star clusters provide natural environments for forming such systems, beyond the limits of standard massive star evolution to core collapse. While repeated BBH mergers can grow BHs through dynamical interactions (the so-called “hierarchical merger” channel), most star clusters with masses ≲10⁶M_⊙have escape speeds too low to retain higher-generation BHs, limiting growth into or beyond the mass gap. In contrast, BH–star collisions with subsequent accretion of the collision debris can grow and retain BHs irrespective of the cluster escape speed. UsingN-body (Cluster Monte Carlo) simulations, we study BH growth and spin evolution through this process, and we find that accretion can drive BH masses up to at least ∼200M_⊙, with spins set by the details of the growth history. BHs up to about 150M_⊙can reach dimensionless spinsχ ≳ 0.7 via single coherent episodes, while more massive BHs form through multiple stochastic accretion events and eventually spin down toχ ≲ 0.4. These BHs later form binaries through dynamical encounters, producing BBH mergers that contribute up to ∼10% of all detectable events, comparable to predictions for the hierarchical channel. However, the two pathways predict distinct signatures: hierarchical mergers yield more unequal mass ratios, whereas accretion-grown BHs preferentially form near-equal-mass binaries. The accretion-driven channel allows dense clusters with low escape speeds, such as globular clusters, to produce highly spinning BBHs with both components in or above the mass gap, providing a natural formation pathway to GW231123-like systems. 

Abstract The dynamical formation of binary black holes (BBHs) in globular clusters (GCs) may contribute significantly to the observed gravitational-wave (GW) merger rate. Furthermore, the Laser Interferometer Space Antenna (LISA) may detect many BBH sources from GCs at mHz frequencies, enabling the characterization of such systems within the Milky Way and nearby Universe. In this work, we use Monte CarloN-body simulations to construct a realistic sample of Galactic clusters, thus estimating the population, detectability, and parameter measurement accuracy of BBHs hosted within them. In particular, we show that the GW signal from 0.7 ± 0.7, 2.0 ± 1.7, 3.6 ± 2.3, and 13.4 ± 4.7 BBHs in Milky Way GCs can exceed the signal-to-noise ratio (SNR) threshold of SNR = 30, 5, 3, and 1 for a 10 yr LISA observation, with ∼50% of detectable sources exhibiting high eccentricities (e ≳ 0.9). Moreover, the Fisher matrix and Bayesian analyses of the GW signals indicate that these systems typically feature highly resolved orbital frequencies (δf_orb/f_orb ∼ 10⁻⁷to 10⁻⁵) and eccentricities (δe/e ∼ 10⁻³to 0.1), as well as a measurable total mass when SNR exceeds ∼20. Notably, we show that high-SNR BBHs can be confidently localized to specific Milky Way GCs with a sky localization accuracy ofδΩ ∼ 1 deg², and we address the large uncertainties in their distance measurement (δR ∼ 0.3–20 kpc). The detection and localization of even a single BBH in a Galactic GC would allow accurate tracking of its long-term orbital evolution, enable a direct test of the role of GCs in BBH formation, and provide a unique probe into the evolutionary history of Galactic clusters. 

Abstract Merging binary black holes (BBHs) formed dynamically in dense star clusters are expected to have uncorrelated spin–orbit orientations since they are assembled through many random interactions. However, measured effective spins in BBHs detected by LIGO/Virgo/KAGRA hint at additional physical processes that may introduce anisotropy. Here we address this question by exploring the impact of stellar collisions and accretion of collision debris on the spin–orbit alignment in merging BBHs formed in dense star clusters. Through hydrodynamic simulations, we study the regime where the disruption of a massive star by a BBH causes the stellar debris to form individual accretion disks bound to each black hole (BH). We show that these disks, which are randomly oriented relative to the binary orbital plane after the initial disruption of the star, can be reoriented by strong tidal torques in the binary near pericenter passages. Following accretion by the BHs on longer timescales, BBHs with small but preferentially positive effective spin parameters (χ_eff≲ 0.2) are formed. Our results indicate that BBH collisions in young massive star clusters could contribute to the observed trend toward small positiveχ_eff, and we suggest that the standard assumption often made that dynamically assembled BBHs should have isotropically distributed BH spins is not always justified. 

Context.The detection of supermassive black holes (SMBHs) in high-redshift luminous quasars may require a phase of rapid accretion, and as a precondition, substantial gas influx toward seed black holes (BHs) from kiloparsec or parsec scales. Our previous research demonstrated the plausibility of such gas supply for BH seeds within star-forming giant molecular clouds (GMCs) with high surface density (∼10⁴ M_⊙ pc⁻²), facilitating “hyper-Eddington” accretion via efficient feeding by dense clumps, which are driven by turbulence and stellar feedback. Aims.This article presents an investigation of the impacts of feedback from accreting BHs on this process, including radiation, mechanical jets, and highly relativistic cosmic rays. Methods.We ran a suite of numerical simulations to explore diverse parameter spaces of BH feedback, including the subgrid accretion model, feedback energy efficiency, mass loading factor, and initial metallicity. Results.Using radiative feedback models inferred from the slim disk, we find that hyper-Eddington accretion is still achievable, yielding BH bolometric luminosities of as high as 10⁴¹ − 10⁴⁴ erg/s, depending on the GMC properties and specific feedback model assumed. We find that the maximum possible mass growth of seed BHs (ΔM^max_BH) is regulated by the momentum-deposition rate from BH feedback,ṗ_feedback/(Ṁ_BHc), which leads to an analytic scaling that agrees well with simulations. This scenario predicts the rapid formation of ∼10⁴M_⊙intermediate-massive BHs (IMBHs) from stellar-mass BHs within ∼1 Myr. Furthermore, we examine the impacts of subgrid accretion models and how BH feedback may influence star formation within these cloud complexes. 

Abstract Dynamical interactions in dense star clusters could significantly influence the properties of black holes, leaving imprints on their gravitational-wave signatures. While previous studies have mostly focused on repeated black hole mergers for spin and mass growth, this work examines the impact of physical collisions and close encounters between black holes and (noncompact) stars. Using Monte CarloN-body models of dense star clusters, we find that a large fraction of black holes retained upon formation undergo collisions with stars. Within our explored cluster models, the proportion of binary black hole mergers affected by stellar collisions ranges from 10%–60%. If all stellar-mass black holes are initially nonspinning, we find that up to 40% of merging binary black holes may have components with dimensionless spin parameterχ ≳ 0.2 because of prior stellar collisions, while typically about 10% have spins nearχ = 0.7 from prior black hole mergers. We demonstrate that young star clusters are especially important environments, as they can produce collisions of black holes with very massive stars, allowing for significant spin-up of the black holes through accretion. Our predictions for black hole spin distributions from these stellar collisions highlight their sensitivity to accretion efficiency, underscoring the need for detailed hydrodynamic calculations to better understand the accretion physics following these interactions. 

Abstract Tidal disruption events (TDEs) are an important way to probe the properties of stellar populations surrounding supermassive black holes. The observed spectra of several TDEs, such as ASASSN-14li, show high nitrogen-to-carbon (N/C) abundance ratios, leading to questions about their progenitors. Disrupting an intermediate- or high-mass star that has undergone CNO processing, increasing the nitrogen in its core, could lead to an enhanced nitrogen TDE. Galactic nuclei present a conducive environment for high-velocity stellar collisions that can lead to high mass loss, stripping the carbon- and hydrogen-rich envelopes of the stars and leaving behind the enhanced nitrogen cores. TDEs of these stripped stars may therefore exhibit even more extreme nitrogen enhancement. Using the smoothed particle hydrodynamics codeStarSmasher, we provide a parameter space study of high-velocity stellar collisions involving intermediate-mass stars, analyzing the composition of the collision products. We conclude that high-velocity stellar collisions can form products that have abundance ratios similar to those observed in the motivating TDEs. Furthermore, we show that stars which have not experienced high CNO processing can yield low-mass collision products that retain even higher N/C abundance ratios. We analytically estimate the mass fallback for a typical TDE of several collision products to demonstrate consistency between our models and TDE observations. Lastly, we discuss how the extended collision products, with high central to average density ratios, can be related to repeated partial TDEs like ASASSN-14ko and G objects in the Galactic center. 

Abstract The existence of compact stellar remnants in the mass range 2–5M_⊙has long been debated. This so-called lower-mass gap (LMG) was initially suggested by the lack of low-mass X-ray binary observations with accretors about 2–5M_⊙, but it has recently been called into question following newer observations, including an LMG candidate with a millisecond pulsar (MSP) companion in the dense globular cluster NGC 1851. Here, we model NGC 1851 with a grid of similar dense star clusters utilizing the state-of-the-art Monte CarloN-body code Cluster Monte Carlo, and we specifically study the formation of LMG black holes (BHs). We demonstrate that both massive star evolution and dynamical interactions can contribute to forming LMG BHs. In general, the collapse of massive remnants formed through mergers of neutron stars (NSs) or massive white dwarfs produces the largest number of LMG BHs among all formation channels. However, in more massive clusters, supernova core collapse can contribute comparable numbers. Our NGC 1851-like models can reproduce MSP—LMG BH binaries similar to the observed system. Additionally, the LMG BHs can also become components of dynamically assembled binaries, and some will be in merging BH–NS systems similar to the recently detected gravitational wave source GW230529. However, the corresponding merger rate is probably ≲1 Gpc⁻³yr⁻¹. 

Abstract Gravitational-wave observations provide the unique opportunity of studying black hole formation channels and histories—but only if we can identify their origin. One such formation mechanism is the dynamical synthesis of black hole binaries in dense stellar systems. Given the expected isotropic distribution of component spins of binary black holes in gas-free dynamical environments, the presence of antialigned or in-plane spins with respect to the orbital angular momentum is considered a tell-tale sign of a merger’s dynamical origin. Even in the scenario where birth spins of black holes are low, hierarchical mergers attain large component spins due to the orbital angular momentum of the prior merger. However, measuring such spin configurations is difficult. Here, we quantify the efficacy of the spin parameters encoding aligned-spin (χ_eff) and in-plane spin (χ_p) at classifying such hierarchical systems. Using Monte Carlo cluster simulations to generate a realistic distribution of hierarchical merger parameters from globular clusters, we can infer mergers’χ_effandχ_p. The cluster populations are simulated using Advanced LIGO-Virgo sensitivity during the detector network’s third observing period and projections for design sensitivity. Using a “likelihood-ratio”-based statistic, we find that ∼2% of the recovered population by the current gravitational-wave detector network has a statistically significantχ_pmeasurement, whereas noχ_effmeasurement was capable of confidently determining a system to be antialigned with the orbital angular momentum at current detector sensitivities. These results indicate that measuring spin-precession throughχ_pis a more detectable signature of hierarchical mergers and dynamical formation than antialigned spins. 

Abstract The origins and mergers of supermassive black holes (SMBHs) remain a mystery. We describe a scenario from a novel multiphysics simulation featuring rapid (≲1 Myr) hyper-Eddington gas capture by a ∼1000M_⊙“seed” black hole (BH) up to supermassive (≳10⁶M_⊙) masses in a massive, dense molecular cloud complex typical of high-redshift starbursts. Due to the high cloud density, stellar feedback is inefficient, and most of the gas turns into stars in star clusters that rapidly merge hierarchically, creating deep potential wells. Relatively low-mass BH seeds at random positions can be “captured” by merging subclusters and migrate to the center in ∼1 freefall time (vastly faster than dynamical friction). This also efficiently produces a paired BH binary with ∼0.1 pc separation. The centrally concentrated stellar density profile (akin to a “protobulge”) allows the cluster as a whole to capture and retain gas and build up a large (parsec-scale) circumbinary accretion disk with gas coherently funneled to the central BH (even when the BH radius of influence is small). The disk is “hypermagnetized” and “flux-frozen”: dominated by a toroidal magnetic field with plasmaβ∼ 10⁻³, with the fields amplified by flux-freezing. This drives hyper-Eddington inflow rates ≳1M_⊙yr⁻¹, which also drive the two BHs to nearly equal masses. The late-stage system appears remarkably similar to recently observed high-redshift “little red dots.” This scenario can provide an explanation for rapid SMBH formation, growth, and mergers in high-redshift galaxies. 

Abstract Alongside the population of several hundred radio millisecond pulsars currently known in Milky Way globular clusters, a subset of six slowly spinning pulsars (spin periods 0.3–4 s) are also observed. With inferred magnetic fields ​​​​​​≳10¹¹G and characteristic ages ≲​​​​​​10⁸yr, explaining the formation of these apparently young pulsars in old stellar populations poses a major challenge. One popular explanation is that these objects are not actually young but instead have been partially spun up via accretion from a binary companion. In this scenario, accretion in a typical low-mass X-ray binary (LMXB) is interrupted by a dynamical encounter with a neighboring object in the cluster. Instead of complete spin-up to millisecond spin periods, the accretion is halted prematurely, leaving behind a “partially recycled” neutron star. In this Letter, we use a combination of analytic arguments motivated by LMXB evolution andN-body simulations to show that this partial recycling mechanism is not viable. Realistic globular clusters are not sufficiently dense to interrupt mass transfer on the short timescales required to achieve such slow spin periods. We argue that collapse of massive white dwarfs and/or neutron star collisions are more promising ways to form slow pulsars in old globular clusters.